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Abstract. We consider a 3rd-order generalized Monge- Ampere equa- 
tion Uyyy — u xxy + u xxx u xyy = (which is closely related to the asso- 
ciativity equation in the 2-d topological field theory) and describe all 
integrable structures related to it (i.e., Hamiltonian, symplectic, and re- 
cursion operators) . Infinite hierarchies of symmetries and conservation 
laws are constructed as well. 
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Introduction 



Monge- Ampere equations [13] is one of the most interesting objects to ap- 
ply methods of geometrical theory of differential equation. Generalizations 
C/3 ! of classical Monge- Ampere equations are discussed, e.g., in [2]. One of such 

£H ■ generalizations is the equation 

U yyy ~ U xxy + u xxxU X yy = 0. (1) 

This is a third-order Monge- Ampere equation ([2j [TO]), but this does not 
{NJ | help too much in understanding its integrability properties. 

Equation (TTJ) is closely related to the associativity equation in 2-d topolog- 
ical field theory |4j and was studied in a number of papers ([5l|6l|8l[9l[II]) and 
' its integrability (existence of a bi-Hamiltonian structure) was established. 

Note though that in these papers the equation was not considered in the 
t^J- | initial form ([I]) , but was rewritten as a three-component system 

a y = b x , by = c x , c y = (b 2 - ac) x (2) 

of hydrodynamical type. Of course, equations ([1]) and ([2]) are closely related, 
but not the same and even not equivalent being associated to each other by 
^ | the differential substitution 

b ■ h - 

1 " — "ran — U XX y: C — ^xyy 

(just like the KdV and mKdV are related by the Miura map or the Burgers 
and heat equations by the Cole-Hopf transformation). 

The aim of this paper is to attack Equation ([T|) directly, not reducing 
it to the evolutionary form, and to study the structures that arise on this 
equation. To this end, we use geometrical and cohomological methods de- 
scribed initially in [TO] and discussed in detail in a recent review paper |12j . 
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These methods has been successfully applied to a number of equation (see, 

e.g., mm. 

In Section[T]we briefly recall basic notions from the geometry of jet spaces. 
Section [2] contains main results on the Monge- Ampere equation (pQ) (includ- 
ing description of Hamiltonian, symplectic, and recursion operators, as well 
hierarchies of symmetries and conservation laws). In particular, we show 
that Equation ([1]) admits a symplectic structure of the form D x (this is the 
only local operator that is responsible for the integrability of the equation 
and it corresponds to the symplectic structure described in [3j [18] ) . A non- 
local Hamiltonian structure D~ l corresponds to this operator. Other oper- 
ators are quite complicated and are described in Sections I2.4.1| 12.4.21 12.5.11 
and I2TTT21 

All computations were done using CDIFF, a REDUCE package for computa- 
tions in geometry of differential equations (see http://gdeq.org). 



1. Theoretical background 

1.1. Jets and equations. Recall that geometric approach to PDEs pQ as- 
sumes that an equation <§ together with all its prolongations (i.e., differential 
consequences) is a submanifold in the manifold J°°(ir) of infinite jets of some 
bundle ir: E — > M, where M and E are smooth manifolds of dimensions n 
and n + m, respectively. 

The first manifold is the one that contains independent variables, while 
the sections of n play the role of unknown functions (fields) in <§ . If U C M 
is a coordinate neighborhood such that 7r|^ is trivial then we choose local 
coordinates x , . . . ,x n in U and u , ... ,u m in the fiber of if\u- Then the 
corresponding adapted coordinates vP a , a being a multi-index, in J°°(7r) are 
defined as follows. For a local section / = (f 1 , . . . , f m ) we set 

dkl p 



dx a 



Functions on J°°(ir) may depend on x % and finite number of u a only. 
The vector fields 

^ d sr^ i 9 . 

JJj = + > u - — - i = l, . . . , n 

are called total derivatives and differential operators in total derivatives are 
called ^-differential operators. 

If an equation is given by the system F = 0, where F = (F , . . . , F r ) is a 
vector-function on J°°(7r), then its infinite prolongation S is given by 

D a {F) = 0, H > 0, 

where D a = D ai o • • • o D (Ts for a = o\ . . .a s . Total derivatives can be 
restricted to $ (we preserve the same notation for these restrictions) and 
generate the Cartan distribution ^ . This distribution is integrable in a 
formal sense, i.e., [X, Y] G ^ for any X, Y € ff, and its n-dimensional 
integral manifolds are solutions of $ . 



INTEGRABLE STRUCTURES FOR A MONGE-AMPERE EQUATION 



3 



1.2. Symmetries. Denote by tToq: £ —¥ M the natural projection. A ir^- 
vertical vector field X on £ is called a symmetry of £ if it preserves the 
Cartan distribution, i.e., if [X,^] C *rf. Every symmetry of £ is of the form 

where summation is taken over internal coordinates on <f and the vector- 
function (p = (ip 1 , . . . , (p m ) satisfies the equation 

= °- 

Here ig is the linearization operator of the vector function F restricted to £ . 
To be more precise, we take the functions F a that define the equation £ 
and construct the matrix ^-differential operator 



i F 



2->cr du\ ■ ■ ■ t^o du™ 

K l^a 9«l G ■ ■ ■ l^a du™ / 



Since thus defined operator is a ^-differential operator, it can be restricted 
to £ and we set 

lg = l F \g. 

The function (p is called the generating function (or section, or character- 
istic) of the corresponding symmetry 9 V and we usually make no distinction 
between symmetries and their generating functions. The set of symmetries 
is a Lie algebra over M with respect to the commutator. We denote this alge- 
bra by sym(<f ). The bracket of vector fields induces a bracket of generating 
functions by 

which is called the Jacobi bracket and is presented by 

W2 H = v:(zM rf )M_ „ W) M) 

in coordinates. 

1.3. Conservation laws. Consider the space A 1 (£) of differential 1-forms 
on £ . It consists of finite sums 



Ai and BJ being smooth functions on £ . The space A 1 (<f) splits naturally 
into the direct sum 

A 1 (£)=A 1 h (£)®A 1 v (£), (3) 

where 

k\{£) = {uj£ A\£) \u = Y,A i dx i } 

is the subspace of horizontal forms, while A\{£) is generated by the differen- 
tial forms ojI = dui — ^ u J ai dx l and is the subspace of vertical (or Cartan) 
forms. 



4 P. KERSTEN, I. KRASIL'SHCHIK, A. VERBOVETSKY, AND R. VITOLO 



Splitting (|3|) generates the splitting 

p+q=s 

where 

A* (#) = A* (*) A ••• A A\ (<?), A*(<?) = A}(0 A ••• A A}(<f ) . 

S v ' * v ' 

g times p times 

Let us introduce the notation A p '' J ((f) = A^(<?) ® A?((f). Consequently, 
the de Rham differential d: A s ((f) — > A s+1 (<f) splits into the sum of the 
horizontal 

4: A p '"(<f ) -> A p ' 9+1 (<^) 

and vertical 

d v : K p ' q {S) -)■ A p+1 ' 9 (<f ) 

parts and one has 

[d fc ,d„] = 0. (4) 

Due to ©, we have a bi-complex structure on A* ((f) which is a particular 
case of Vinogradov's "^-spectral sequence [151 EH fT7] . Denote by E\ ,q {S) 
the cohomology of dh at the term AP ,q {S ). Then d v induces the differentials 

The group E®' 71-1 ^) plays a special role in the theory. Its elements are 
called conservation laws of $ , while the group itself is denoted by Cl(S'). 
We also shall need the group £'J' n " 1 whose elements are called co symmetries 
and which is denoted by cosym(<f). 

To proceed, we shall need additional constructions. Let P and Q be the 
spaces of sections of vector bundles over S. Let P = Hom(P, h n {<§)) and 
similar for Q. Then for any "^-differential operator A: P — > Q its formally 
adjoint A* : Q — > P is defined by the Green formula 

(A*(q),p) - (q,A(p)) = d h uj(q,p), 

where u: Q x P — >• A™" 1 ((f) is a map which is a "^-differential operator in 
both arguments and (• , •) denotes the natural pairing. If A is given by the 
matrix (Ay), where Ay = J2a a ij-Dcr, then A* = (A*^), where 

In what follows, we shall assume that equation at hand satisfies the fol- 
lowing conditions: 

(1) The differentials dF 3 of the functions that define $ are linear inde- 
pendent at all points of $ . 

(2) If A is a ^-differential operator such that A o tg = then A = 0. 

(3) If A is a "^-differential operator such that A o £J = then A = 0. 

If an equation enjoys these conditions then the following statements are 
valicS 



They follow from Vinogradov's 2-Line Theorem [17| . 
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(1) The differential 5: Cl(£) — > cosym^) is monomorphic, i.e., S(lo) = 
if and only if oj = 0. 

(2) The group of cosymmetries coincides with the kernel of i.e., ip € 
cosym((a') if and only if £g(ip) = 0. 

If u> £ Cl(<^) is a conservation law then the cosymmetry S(lo) is called its 
generating function (or generating section). 

1.4. Differential coverings. Let £ and £ be two equations. A smooth 
map r : — > £ is called a morphism if it takes the Cartan distribution on <f 
to that on £ '. A surjective morphism t is said to be a covering if for any 
point 6 € £ the differential dr| e maps the Cartan plane < tfe(£) to ^ T ^(£) 
isomorphically. Coordinates along the fibers of r are called nonlocal variables 
in the covering under consideration. Let r' : £' — > £ be another covering. 
We say that it is equivalent to r if there exists a morphism / : £ ^ £' which 
is a diffeomorphism and such that r = t' o /. 

If Z?i, . . . , D„ are total derivatives on and u; , . . . , w r , . . . are nonlocal 
variables then the covering structure is given by vector fields 

Di = Di + X i} i = l,...,n, (5) 

where Xi = ^2 a Xfd/dw a are r- vertical fields that satisfy the condition 

Di(Xj) - D 3 (Xi) + [X h Xj] = 0, 1 < i < j < n. (6) 

A covering is Abelian if the coefficients Xf do not depend on nonlocal vari- 
ables. In this case, d6]) amounts to 

Di(Xj) - Dj(Xi) = 0, l<i<j<n. (7) 

In the particular case of one-dimensional coverings, conditions ([7]) define 
a dh-closed horizontal 1-form uj t = ]T\ Xi dx l on £ and two coverings of 
this type are equivalent if and only if the corresponding forms are in the 
same cohomology class, i.e., u T — uj t i = dh(g) for some function g. When n 
(the number of independent variables) equals two, this establishes a one-to- 
one correspondence between the group Cl(<#) and the equivalence classes of 
one-dimensional Abelian coverings over £ . 

If t: $ — >• $ is a covering then symmetries of <§ are called nonlocal r- 
symmetries of £ '. Note also that any ^-differential operator A on £ can be 
lifted to a ^-differential operator A on £ . This is being done by changing 
total derivatives Di in the local representation of A to Di using ([5]). In 
particular, the linearization operator ig can be lifted in such a way and 
solutions of the equation 

e*(<p) = 

are called (nonlocal) shadows of symmetries in the covering r. In a similar 
way, solutions of 

~W) = o 

are (nonlocal) shadows of cosymmetries in the covering r. 



6 P. KERSTEN, I. KRASIL'SHCHIK, A. VERBOVETSKY, AND R. VITOLO 

1.5. The ^-covering. Let £ be an equation. Consider a new set of de- 
pendent variables q = (q 1 ,...,q m ) (in many respects it is convenient to 
consider q as an odd variable), where m is the number of unknown functions 
in £ , and augment the initial equation with 

ts(q) = 0. (8) 

The resulting system, consisting of £ and equation (JSD, is called the l- 
covering of £ '. It is an analog of the tangent bundle for the equation £ . The 
^-covering is important for the subsequent computations due to the following 
properties. 

1.5.1. Recursion operators for symmetries. Consider a vector-function <3? = 
($ , . . . , $ m ), where = <&a i0 -(/", which is a symmetry shadow in the 
^-covering. This means that it satisfies the equation 

M$) = o, (9) 

where lg is the linearization operator lifted to the ^-covering. Then it can 
be shown that the matrix ^-differential operator ^$ = (%2 ff &a,aD a ) takes 
symmetries of £ to symmetries. In other words, ^$ is a recursion operator 
for symmetries. 

A recursion operator S% is called hereditary if 

{Mifx ,Mip 2 }- M{{Mip x , <p 2 } + { <pi , &<p 2 } - ><%{V2 , }) = 

Hereditary operators possess the following property important for integra- 
bility: let (p be a symmetry such that 

9 v (£)-[£ v ,O\=0. 

Then all symmetries M l (p pair-wise commute, i.e., form a commutative hi- 
erarchy. 

1.5.2. Symplectic operators. In a similar way, let us consider now a vector- 
function ^ = (ty 1 , . . . where = Eaff*«. ff 9" (recall that r is the 
number of the functions F 3 that define £) , that satisfies the equation 

= o, (10) 
where £J is the lift of the operator l* s to the ^-covering. Then the opera- 
tor S% = (Yla ^a,aDa) takes symmetries of the equation 3 to its cosym- 
metries. 

Let an operator 5? satisfy the condition 

y*olg = t g oy. (11) 

Then 5? is identified with a variational 2- form Qy on £ whose values on 
symmetries is given by 

This form can be considered as an element of the group Ef' n ~\£) in the 
term E\ of the ^-spectral sequence. 

Let now u\, 002 be two conservation laws such that 

5uJi = J?<pi, i = 1,2, 
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for some (pi, (f2 G sym<<?. Then the bracket 

{u)x,u) 2 }y = n^(yi,y2) 

is defined. This bracket is skew-symmetric by (jlip and satisfies the Jacobi 
identity if 

sn y = o. (12) 

Operators that enjoy properties (fTT]) and (fT2|) are called symplectic. 

1.5.3. Nonlocal covectors. Solving equations ([9]) or (fTD]) leads often to trivial 
results only. The reason is that symplectic and especially recursion operators 
are in many cases nonlocal, i.e., contain terms like D^ 1 . Such terms are 
incorporated into solution by introducing nonlocal variables that amounts 
to constructing appropriate coverings. One of the ways to construct the 
latter is based on the following fact (see [12]): to any cosymmetry of the 
equation £ there corresponds a conservation law on the l-covering. We 
call these conservation laws nonlocal covectors. Consequently, if n = 2 an 
Abelian covering corresponds to a cosymmetry. Numerous computations [7J 
\10\ fTT] show that nonlocal variables arising in such a way are sufficient to 
find necessary structures. 

1.6. The ^-covering. Let again £ be an equation. Consider a another 
new set of dependent variables p = (p 1 , . . . ,p r ), where r is the number 
of functions F J that determine the equation £, and augment the initial 
equation with 

t g {p) = 0. (13) 

The resulting system, consisting of £ and equation (fT3j) . is called the l-- 
covering of £ . It is an analog of the cotangent bundle for the equation £ . 
The ^""-covering is also important for the subsequent computations due to 
the following properties. Like the variable q in the ^-covering, it is convenient 
to consider the variable p to be odd. 

1.6.1. Hamiltonian operators. Consider a vector-function = (Q 1 , . . . , $ m ), 
where & = $a i(7 p", and assume that it is a symmetry shadow in the 
^""-covering. This means that it satisfies the equation 

M$) = 0, (14) 

where lg is the linearization operator lifted to the ^-covering. Then it can 
be shown that the the matrix ^-differential operator = (Yla ^ce,aD a ) 
takes cosymmetries of £ to symmetries. 

Solutions to (|14p of special type are identified with variational bivec- 
tors Ajf on £ . These solutions must satisfy the condition 

lgoM? = ,3^* otg. 

In this case, the operation 

{wi,w 2 }jr = (Jf(6(jJi),5u)2), ui,u> 2 G Cl(£), 

defines a skew-symmetric bracket on the space of conservation laws. This 
bracket satisfies the Jacobi identity if and only if {Ajg>, A,^»] = 0, where [•, •] 
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is the variational Schouten bracket (see [10\ I12j) on the space of variational 
multi- vectors. In this case, one has 

1.6.2. Recursion operators for co symmetries. Finally, let us now consider a 
vector-function ^ = (vp , . . . , ^ r ), where flfi = ^ Q ^Q )0 -p°, and assume 
that it is a cosymmetry shadow in the ^-covering. This means that it 
satisfies the equation 

= 0, (15) 
where lg is the linearization operator lifted to the £*-covering. Then it can 
be shown that the the matrix ^-differential operator My = (^ CT ^a^D^) 
takes cosymmetries of $ to cosymmetries. In other words, ffl is a recursion 
operator for cosymmetries of £ '. 

1.6.3. Nonlocal vectors. Similar to Section [T31 equations (fbi]) and (fT5|) lead 
often to trivial results only. The reason is the same: Hamiltonian and re- 
cursion operators are in many cases nonlocal. Such terms are incorporated 
into solution by introducing nonlocal variables that amounts to constructing 
appropriate coverings. One of the ways to construct the latter is based on 
the following fact: to any symmetry of the equation £ there corresponds a 
conservation law on the £* -covering. We call these conservation laws non- 
local vectors. Consequently, if n = 2 an Abelian covering corresponds to a 
symmetry. As computations [TJ [TUT, [TT] show, nonlocal variables arising in 
such a way are sufficient to find necessary structures. 

1.7. General computational scheme. In all the computations we did 
to analyze particular equations (Equation (TTJ included) we adhered to the 
following scheme: 

• Extension of the initial equation with a minimal set of nonlocal 
variables (usually associated with conservation laws) to ensure exis- 
tence of nontrivial solutions to the main equations defining integrable 
structures. 

• Computation of a minimal set of (local and nonlocal) symmetries 
and cosymmetries necessary to (a) hierarchy generation and (b) con- 
struction of nonlocal vectors and covectors. 

• Extension of the ^-covering and and construction of symplectic struc- 
tures and recursion operators for symmetries. 

• Extension of the ^-covering and and construction of Hamiltonian 
structures and recursion operators for cosymmetries. 

2. Main results 

For internal coordinates on $ we choose the functions 

d k+i u 

Uk,i = dx kQ y i ' i = 0,1,2, k = Q, 1,..., (16) 
and then the total derivatives on $ acquire the form 

D x = — + y7ufc+i,o q 1" Uk+1,1 q h Uft+1,2-5 ), (17) 

ox ■ du kfi du kl du k2 
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fc>0 

Equation (pQ) is homogeneous with respect to the following weights 
\u\ = 0, \x\ = —1, \y\ = —4. 

2.1. Conservation laws and Abelian coverings. In the sequel, we shall 
need nonlocal variables that will be denoted by Qij. The second subscript 
here indicates the weight of the variable, while first the one corresponds to 
the level of nonlocality. By the latter we mean the following. The variables 
of zero level are determined by local functions on <§: 

dQoj . dQ ,7 

= -^0,1^4,0 + 1*0,2, — a = - U 0,1U 3 1 - Mo,2^3,0 + 1*1,11*2,1, 

ox oy 

9Qo,9 ,2 2 9Q g 2 

—q— = 21*0,2^2,0 + - ^2,0^2,1, q = 2no,2«l,l - ^1,2«2,0) 

and 

dQo,i2 , s <9<2o,12 1/ 2 o i 

— ^— = ^^^0,2 - lt2,0«2,l), q = 2l U 0,2 - 2,1ii i i,Ul j 2«2,o)- 

The variables of level 1 are determined by local functions and by variables 
of zero level: 

9 Ql,3 2 °Qi,3 ofri . x 

^ = 2li ,l - « 2) 0) = 2{Qo,7 + «0,l'"3,0 - Wl,l«2,0j) 

9Qi,e „ , <9Qi,6 1. 2 



and 



ggij 
<9y 



Qo,7 + 1*0,1«3,0, a ' = o U l 1' 
Ox Oy 2 ' 



: Qo,9 - 2-uo,ini i0 n 4i o - 4u ,ii*20i*30 + 2uo,2«l,o> 

: 4Q ,12 - 2«0,l«l,0'"3,l - 2lt ,l«l,l«3,0 - 2uo,lM2,0^2,l 
2«0,2l*l,0l*3,0 - 1*0,21*2,0 + 2Ul,o1*l,l«2,l + 2l/ 1 , 1 ti2,0- 



There exist deeper level nonlocalities, such as 

- -18Qi, 6 - 2u ,in 2i o + 4«i, i*i,i + u 2 , 



dx 
dQ 2 , 5 
dy 



-3Q ,9 - 2ii ,ii*i,i + 4tt ,2i*i 







and 



0Q2 7 2 
q x = -40Qo,7«i i0 - 10Qi, 8 - 10w 0>1 - 60n ,ini i0 ii3,o 

+ «1 ) 0«2,0 U 3,0 - 7}U2 , 
0Q2 7 2 

- -40Q ,7i*oi - 10Qi,ii ~ 30u i« 3) o + 20uo,i«i,ii*2,o 



+ 10tt ,11i| 1*3,0 - 30ltl )O lti I + 1tl,0l*2,0 n 2,l - 3«l,l1i| o, 
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as well as 

— « — = 0^2,5 + Qi,3U 2 ,o - -n ,iMi,o, 
ox 3 3 

d( ^y A = 2 Qoj u i,o + Ql,3^1,l - Ql,8 - 2uo,i ~ u 0,lU 2 fl- 

Remark 1. The zero level nonlocal variables are associated to conservation 
laws of equation (P). For example, to Qo,7 there corresponds the conservation 
law 

^0,7 = (-^0,1^4,0 + "0,2) dx + {-u 0> iu^i - u ,2M3,o + ui t xu 2 ,i) dy. 

The first level nonlocal variables are associated to conservation laws of the 
Abelian coverings determined by the zero level variables, etc. 

Remark 2. Of course, the list of nonlocal variables above is not exhaustive 
at all. We described only those ones that are used to construct the necessary 
nonlocal symmetries (§[22} and cosymmetries (§ I2.3p ; see also § 11.71 New 
nonlocalities do arise under the actions of recursion operators, but there is 
no need to describe them explicitly here. 

2.2. Symmetries. This direct computation is needed by the two reasons: 
to construct nonlocal vectors (see Subsection II .6.3| ) and to use the obtained 
symmetries as 'seeds' for the hierarchies. 

The linearization of equation fl]) has the form 

Dy(<P) ~ 2u 2il D 2 x D y (<p) + u lj2 D 3 x (<p) + u 3fi D x D 2 y &) = 0, (18) 

where the total derivatives D x and D y are given by (|17p . 
Solving fig]) , we get the following solutions^. 

2.2.1. Symmetries of degree 0: 

<P0 = 1, ipi=Ui t0 , ipl = U ,l, = Q2,5 + 8u ,lUi ;0 , 

ipl = -2Q ,7«i,o + Qi,8 - 2no 5 i + «o,i'"|,o- 

2.2.2. Symmetries of degree 1: 

V-4 = V, P-i = x, 
¥>o,i = xui fi - 4u, <pl,2 = y u o,i + u > 
tp\ = 4xu ,i - Ql,3, 
(f\ = x(Q 2 ,5 + 8«o,i u i,o) - 3(8(5 3>4 - 3Qi j3 tti )0 + 16tttt 0) i). 

2.2.3. Symmetries of degree 2: 

ip\ = x 2 u lfi + 4xyn ,i - Axu - yQi, 3 , 

2 _ o„,2„._ . ,.2 



Lp 2 = 2x Uq,1 - xQi :3 - u 



i,o- 



9 

In the notation for symmetries, the superscript indicates the polynomial order of a 
symmetry with respect to x and y, the first subscript equals the weight, while the second 
one, if any, is the number of a symmetry in the set given by particular weight and order. 
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2.2.4. Symmetries of degree 3: 

<P-3 =x 3 - 2yu lfi , 
<Pl = 12x 3 n ,i - 9x 2 Qi y3 - 18xul - 2y(Q 2 ^ + 8uo,it*i,o) + 24m*i )0 . 

2.2.5. Symmetries of degree 4: We shall also need one symmetry of order 4, 
which is of the form 

f-4 = x 4 — Sxyuifl — 8y 2 iio,i + 16yu. 

2.3. Cosymmetries. The reasons to compute cosymmetries explicitly are 
similar to those indicated in Subsection 12.21 

To find cosymmetries, we are to solve the equation adjoint to (fTSjl. i.e., 

Djty) - 2D 2 x D y (u 2 ^) + D 3 x ( Ul ^) + D x D 2 y (u 3 ^) = 0. (19) 

Using the notation similar to the one from Section 12.21 let us some compu- 
tational results needed below. 

2.3.1. Cosymmetries of degree 0: 

^0 = h ^2= u 2,0, ^5=^1,1) 

^6 = -18Qi, 6 + 6n 0i in 2i o + 12ui i0 iii,i + ti| , 
^9 = 2Q ju 2 ,o ~ Qo,9 + 41*0,1^1 + 2-uo,in 2i on 3i o - wi,i«2,o- 

2.3.2. Cosymmetries of degree 1: 

lp-4 = V, fp-i = x, Vi,i = xu 2fi - 3«i, , 
ipl )2 = 2/1*1,1 + ui fi , tp\ = 4xni ; i + 2u ,i + u 2 2fi . 

2.3.3. Cosymmetries of degree 2: 

ip-2 = 3x 2 - 2yu 2fi , 
ipl = x 2 u 2)0 + ixyu^i - 2xui fi + y(2i*o,i + «i,o) _ 4-u > 
itfl = 2x 2 ui ) i + x(2« ,i + «i,o) ~ <9i,3 - 2u 1)0 u 2fi . 

2.3.4. Cosymmetries of degree 3: 

-3 



= 2m ny - 2uiiy 2 - 2u 2fl xy + x 3 , 



(3 3 \ /3 

-ito.i + ^u| j - x ( -Qi )3 + 3ui i0 u 2 ,o 

+y ^3Qi j6 - u ,ii*2,o - 2i*i,oiti,i - ^2,0 J + (2uu 2fi + ^uf fl ). 

2.4. The ^-covering. This covering is determined by the system of equa- 
tions 

^"yyy ^xxy ^xxx^xyy — 0, 

Qyyy 2u X xyQxxy ~i~ U X yyQxxx ^xxxQxyy — 0, 

where q is an odd variable along the fiber of the covering. 

Internal coordinates on the space of the ^-covering are functions f|16j) 
together with 

Qk+i 

qk ' i= dx^' * = M,2, fc = 0,l,..., (20) 
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while the total derivatives in these coordinates take the form 
d \ - / d d d 



dx 7~^n ^ ' du k,o ' duk,i ' du k 2 

k>0 

d d d 



+ Qk+1,0 q 1" Qk+l,l q 1" Qk+1,2 



OQkfl OQk,i OQk, 



^9v-^/ d d k 2 .9 



^ ^ ^ ' 5 «fc,0 ' ^Mjfe i X ' ' ' du k>2 

k>0 

+ 9fc,l p~~ 1" ?MF 1" £>z(2u2,l<?2,l - "1,293,0 - «3,0?1,2)- 



<-><//,■.(> ' ''' "" ' " " ' ' <><lk.2- 

By the general theory [12], to any cosymmetry ip of the initial equation 
there corresponds a conservation law on the ^-covering (a nonlocal form). 
Denote by f2w, the corresponding nonlocal variable. In the case of Equa- 
tion (fTjh this variable is defined by the relations 



dx 
dy 

where 



and 



V"?o,2 + ao,i<?o,i + aofiQ, 

^0,250,2 + + ^2,092,0 + &o,i?o,l + &l,o<?l,o + ^o,o9, (21) 

bo,2 = -u 3fi ip, &i,i = 2n 2 ,iV', &2,o = -ui,2ip, (22) 
6 ,i = --Dx(&i,i) 5 fo i,o = --Dx(&2,o), 
^o,o = -D x (b 1)0 ) 

«o,i = D x (b ,2) - D y (ip), a ,o = ^(bo.i) - Dj,(oo,i). (23) 
Below we use the notation fl k ■ = fl , h 

2.4.1. Recursion operators for symmetries. To find recursion operators for 
symmetries of Equation ([1]) , we solve the equation 

= 0, 

where lg is the linearization operator (|18p in the ^-covering extended by 
nonlocal forms and $ is a function on this extension. 
The simplest nontrivial solution is of the form 

= n% - x n 2 _ 2 + 4yii\ j2 + 2 y n\ 1 + ^si 1 ^ 

- 2u lj0 rL-4 - 2y 2 n° 5 - 2xySl% + {2u 1>0 y - x 3 )Q° . (24) 

Using the first of Equations (|2ip . we put into correspondence to every 
nonlocal form Q,^ the operator 

% = D' 1 o (^£>J + oo,i D y + oo,o) (25) 
where the coefficients are determined using relations (|22p and (j23|) . i.e., 

a ,i = -D x (u 3)0 ip) - D v (ip), 
a , = -2Dl(u 2 ,iip) + D x D y (u 3fi 4>) + D^ip). 
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Then the recursion operator 

#1 = %s_ 3 - x%2_ 2 + 4y^i 2 + 2y^i j + 3x 2 ^ Li 

- 2m i0 ^i 4 - 2y 2 ^o - 2xy^o + (2wi )0 y - x 3 )^o (26) 
corresponds to solution ([24 



Remark 3. Since the variables x and y in Equation (pQ) "enjoy equal rights", 
one can put into correspondence to the nonlocal form fl^, using the second 
equality in (f2Tj) . the operator 

% = D x l o {b^ 2 D 2 y + b 1A D x D y + b 2fi D 2 x + b 0A D y + b lfi D x + b ,o) 

and construct a recursion operator similar to operator ()26p . 

The action of operators £%\ and ffl'^ on symmetries of Equation ([1]) is the 
same. 

2.4.2. Symplectic structures. To find symplectic structures, we solve the 
equation 

= o, 

where, similar to Section 12.4.11 l*g is the operator (|19p on the ^-covering 
extended by the nonlocal forms, while \& is a function on this extension. 

Below we present the first two solutions of this equation. The simplest 
one is 

*i = ni )0j 

and the corresponding symplectic structure 

y 1= D x . (27) 

The next solution is nonlocal: 

$ 2 = n 2 _ 1 - ftcnLi + 2u2, n 1 _ i + 2 y nl - (2u 2 ,oy - 3x 2 )n%. 

The corresponding symplectic operator S? 2 : sym $ —> cosym $ is 

y 2 = - 6x^ Li + 2u 2 ,o%i 4 + 2 y%° 2 ~ ( 2n 2,oy - 3x 2 )%o, 
where the operators 5fy are defined by Equation (|25j) . 

2.5. The £*-covering. As it was indicated above, this covering is obtained 
by adding to the initial equation another one, which is adjoint to the lin- 
earization. In other words, this covering is described by the equations 

^"yyy ^xxy ^xxx^xyy — 0, 
UxxyyPxx 2u X xxyPxy ~i~ ^xxxxPyy 

'UxyyPxxx 2u X xyPxxy ~i~ ^xxxPxyy ~i~ Pyyy — 0] 

where p is a new odd variable. 

Internal coordinates in the space of the -P-covering are functions ()16p 
together with the functions 

Pk,i = a k a v ^ = 0,1,2, k = 0, 1, . . . , 
ox K oy l 

while the total derivatives in these coordinates are of the form 
d \ - ( d d d 

L>x = + > A Uk+lfi q h «fe+l,l q 1" q 

^ V du K0 du k,l OU k , 2 

k>0 
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d d d 

+ Pk+1,0^ 1" Pk+1,1 o 1" Pfc+1,2 



OPkfl OPkA OP, 



k.2 



k>0 

9 d T^kr 

+ Pfc,l a \~Pk,2-^ D x (U 2 ,2P2fl ~ 2lt 3 1 + 1i4,0£>0,2 

OPkfi OPk,l 

+ «1,2P3,0 - 2« 2) lP2,l + U 3) oPl,2) 



9pfc,2 ' 

Let y be a symmetry of Equation (pQ). Then (see [12J) a conservation 
law on the £*-covering corresponds to this symmetry and, consequently a 
nonlocal variable, which we denote by 11^ and call a nonlocal vector. For 
Equation ([1]), the correspondence (p i-> 11^, is given by the relations 

-q^- = W,2 + oo,iPo,i + ao,oP, 

dn 

—r^- = b ,2P0,2 + + &2,0P2,0 + &o,iPo,i + &l,OPl,0 + & o,oP, (28) 

where 

&0,2 = = 2U2.lV, b 2,0 = -«1,2V, 

6 ,1 = --Da;(fel,l) + 2u 3i iv?, bi t0 = -D x {b 2 fl) - U 2<2 ip, 

b ,o = ~D x (b lfi ) (29) 

and 

«0,i = D x (b 0t2 ) - D y (<p) + u ifi (p, a , = Ac(&o,i) - D y (a ,i). 
We use the notation Tlf- = II * below. 

To describe the subsequent results, let us put into correspondence to a 
nonlocal vector H v the operator 

^ = D~ l o (pflj + a ,i^ y + ao,o), (30) 

see the first of Equations ([28]) . Here, due to relations ([29]) and ([30]) . the 
coefficients ao,o and ao,i are of the form 

a ,o = -2u 2i iL> 2 (</?) + u 3fi D x D y ((p) + £>y(y?) - U3,iAc(<p), 
ao,i = -U3,oD x {ip) - Dy(tp). 

2.5.1. Hamiltonian structures. Similar to Sections 12.4.11 and \2A.2\ Hamil- 
tonian structures are solutions of the equation 

= 0, (31) 

where the operator ig is the linearization lifted to the £*-covering extended 
by nonlocal vectors. 

The simplest solution of Equation (|3"Tj) is 

$ = n^ 8 -2yni 4 + y 2 n{], 

to which the operator : cosym<f — > sym<f 
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corresponds. The next solution is much more complicated and has the form 
$1 = ni 4 - 4xU 3 _ 3 + 6x 2 II 2 2 - 8ui i0 n^ 5 - 8u ,iII 2 8 + 16yII^ 2 
+ (8xu lfi + 16yu ,i - 16u)nl 4 + SyU 1 ^ - (4x 3 - Syui^U 1 ^ 
- 8y 2 n° 4 - 8xyII? + (x 4 - 8xyu lfi - 8y 2 u ,i + 10yu)U% 



and the operator 



3HP X = % U - 4x% 3 _ 3 + &r% a _ a - 8«i j0 V B " 8u oaK-s + 1%^ >2 
+ (8xm i0 + 16yuo,i - 16n)^ L4 + - (4x 3 - 8yui )0 )^i x 

- 8y 2 ^o - 8sy^o + (x 4 - 8xyui, - 8y 2 -u ,i + 16yu)%o. 
corresponds to this solution. Here and below the operators ^ are defined 



by Equality (j30D . 

Remark 4. The form of the above presented solutions is determined by the 
choice of nonlocal vectors in the £*-covering. They, in turn, are due to a 
basis in the space of symmetries. Actually, Equation has a simpler 
solution <3?q = lip where the nonlocal variable n 4 is defined by the system 

^= P ^i = n', 

dx ' dy 2 ' 

<9IT 2 dUL „. 



<9x ' ay 

an 3 an', 

-5— = P0,2, = -^3,0^0,2 + 2« 2 ipi i - Ul 2P2.0- 

ax ay 
To this solution there corresponds the Hamiltonian operator 

which is the inverse to the symplectic operator (|27p . This operator cor- 
responds to the Hamiltonian operator Jo from [S]. The operator M\ that 
explicitly depends on x and y seems to be new. 

Remark 5. There is a relation 

between the two Hamiltonian structures, where is the recursion operator 
given by Equality 



2.5.2. Recursion operators for cosymmetries. To conclude our study, we con- 
sider finally the equation 

= o 

on the ^""-covering extended by nonlocal vectors. Here is its simplest non- 
trivial solution 

* = n 3 3 - 3xn 2 2 + 2u 2j0 n 2 5 + 2u hl n 2 _ s 

+ 2(m, + 2ui.ii/ - « 2 , x)ni 4 - (2u 2 , y - 3x 2 )n^ + 2yU° 1 
- (2u lfi y - 2ui 5 iy 2 - 2u 2 flxy + x 3 )!^. 
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To this solution there corresponds the recursion operator Rq : cosym <f — > 
cosym £ of the form 

Ro = % 3 _ s - 3x% 2 _ 2 + 2«2,o#J_ B + 2u 1A %2_ g 

+ 2(n lj0 - 2u 1>iy - u 2 ,qx)9 v x_ a - (2u 2 , y - 3x 2 )^ + 2y^o 

- (2«i,o2/ - 2u 1>1 y 2 - 2u 2 fixy + x 3 )^o. 

Remark 6. In the case of evolution equations, existence of a commutative 
hierarchy means that the initial equation possesses "higher analogs" that 
are obtained by the action of a recursion operator. In the general case, 
there exists no well defined action of recursion operators on the equation 
and thus such a construction is impossible. Nevertheless, one can consider 
the following alternative scheme: (1) to pass, if possible, to the evolutionary 
presentation of the equation at hand, (2) to construct "higher analogs" in 
this presentation and (3) return back to the initial variables. We did not set 
the question whether such a scheme is invariant (i.e., whether the result is 
completely determined by the initial equation or depends on the procedure). 
Probably, an answer to this question will lead to a deeper understanding of 
integrability for general equations. 
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